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Abstract 

We present a charge and self-energy self-consistent computational scheme for correlated systems 
based on the Korringa-Kohn-Rostoker (KKR) multiple scattering theory with the many-body ef- 
fects described by the means of dynamical mean held theory (DMFT). The corresponding local 
multi-orbital and energy dependent self-energy is included into the set of radial differential equa- 
tions for the single-site wave functions. The KKR Green's function is written in terms of the 
multiple scattering path operator, the later one being evaluated using the single-site solution for 
the t-matrix that in turn is determined by the wave functions. An appealing feature of this ap- 
proach is that it allows to consider local quantum and disorder fluctuations on the same footing. 
Within the Coherent Potential Approximation (CPA) the correlated atoms are placed into a com- 
bined effective medium determined by the dynamical mean field theory (DMFT) self-consistency 
condition. Results of corresponding calculations for pure Fe, Ni and Fe^Mi^ alloys are presented. 



PACS numbers: 71.15.Rf, 71.20.Be, 82.80.Pv, 71.70.Ej 



I. INTRODUCTION 



One of the first band structure methods formulated in terms of Green's functions is the 
KKR method of Koringa, Kohn and Rostoker . Although it is not counted among the 
fastest band structure methods, it is usually regarded as a very accurate technique. The 
advantage of the KKR method lies in the transparent multiple scattering formalism which 
allows to express the Green's function in terms of single-site scattering and geometrical 
or structural quantities. A second outstanding feature of the KKR method is the Dyson 
equation relating the Green's function of a perturbed system with the Green's function of 
the corresponding unperturbed reference system. Because of this property, the KKR Green's 
function method allows to deal with substitutional disorder including both diluted impurities 
and concentrated alloys in the framework of the Coherent Potential Approximation (CPA) 
0| . Within this approach (KKR-CPA) the propagation of an electron in an alloy is regarded 
as a succession of elementary scattering processes due to random atomic scatterers, with an 
average taken over all configurations of the atoms. This problem can be solved assuming 
that a given scattering center is embedded in an effective medium whose choice is open and 
can be made in a self-consistent way. The physical condition corresponding to the CPA is 
simply that a single scatterer embedded in the effective CPA medium should produce no 
further scattering on the average. A similar philosophy is applied also when dealing with 
many-body problems for crystals in the framework of the so called dynamical mean field 
theory (DMFT, for review see Ref. [l). Thus it seems to be rather natural to combine the 
DMFT and KKR methods to arrive at a very reliable and flexible band structure scheme that 
include correlation effects beyond the standard local density (LDA) or generalized (GGA) 
approximations. In fact the combination of the KKR-CPA for disordered alloys and the 
DMFT scheme is based on the same arguments as used by Drchal et al. when combining 
the TB-LMTO Green's function method for alloys [y] with the DMFT. In contrast to their 
approach, however, the formalism presented below allows to incorporate correlation effects 
via a corresponding self-energy when calculating the electronic single-site wave functions. 

First attempts to achieve a self-consistent description of local correlation effects in crys- 
tals have been made already many years ago. In the third paper of a famous series Hub- 
bard |7| has introduced an alloy analogy and by an appropriate decoupling scheme for the 
Green's function a set of equations has been derived that represent a self-consistent formu- 
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lation equivalent to the CPA approximation. In contrast to the DMFT the "Hubbard III" 
approximation considers quantum on-site fluctuations as static ones which leads to some 



iquid properties, missing of the so called 
. Keeping in mind the above conceptual 



shortcomings such as violation of some Fermi 
Kondo peak near the metal-insulator transition 
analogies it is our purpose to present here a combined Local Density Approximation and 
Dynamical Mean Field Theory (LDA+DMFT) electronic structure technique, including the 
case of disordered solids, in the framework of the KKR method. The many-body correlation 
effects are treated by means of the DMFT, while the disorder is described in the frame- 
work of the CPA. Taking into account the local nature of the DMFT approximation the 
self-energy is represented by a local complex energy dependent quantity (which is a matrix 
in orbital indices) viewed as a contribution to the electronic potential. We note that for 
a general non-local energy dependent potential multiple scattering theory offers a solution 
known as the optical potential 8]. However, the nonlocal self-energy is far too complicated 
to be used in a realistic computation. 

Very recently a combined LDA+DMFT computational scheme was proposed in which 
;he so called Exact Muffin-tin band structure method was used. In the EMTO approach 
^, [h| the one-electron effective potential is represented by the optimized overlapping 
muffin-tin potential which is considered as the best possible spherical approximation to the 
full-one electron potential. In essence the one-electron Green's function is evaluated on a 
complex contour similarly to the screened KKR technique, from which it was derived. In 
the iteration procedure the LDA+DMFT Green's function is used to calculate the charge 
and spin densities. Finally, for the charge self-consistent calculation one constructs the new 
LDA effective potential from the spin and charge densities Q], using the Poisson equation 
in the spherical cell approximation Q|. 

In contrast to the EMTO implementation Q], the present work follows a natural de- 
velopment in which the self-energy is added directly to the coupled radial differential equa- 
tions which determine the electronic wave function within a potential well and this way the 
single-site t-matrix. Because this way also the scattering path operator of multiple scatter- 
ing theory used to set up the electronic Green's function is determined unambiguously, no 
further approximations are needed to achieve charge self-consistency. 

The paper is organized as follows. Section |H] presents a general formulation of the prob- 
lem. Section Til Al provides an extension of the derivation of the multiple scattering Green's 
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function to include the self-energy, and in particular provides the information on the many- 
body solver. Section III Bl describes the many-body solver used in our calculation, that is 
based on a modified fluctuating exchange interaction approximation. The combined self- 
consistency cycle is presented in section 111 Ul Finally, results and discussions are presented 
in section UTTl 



II. FORMULATION OF THE PROBLEM 



The DMFT method has already been implemented within several band structure methods 
based on a wave function formalism: first in the linear muffin-tin orbital method in atomic 
sphere approximation (ASA-LMTO) 14, 15, l3| and then in full-potential LMTO 



as well as in a screened KKR or exact muffin-tin orbitals approach (EMTO)[l2j. The 
emerged LDA+DMFT method can be used for calculating the electronic structure for a 
large variety of systems with different strength of the electronic correlations (for a review, 




see Refs. llflL \2(f). To underline the importance of complete LDA+DMFT self-consistency 
we mention that the first successful attempt to combine the DMFT with LDA charge self- 
consistency gave an important i nsig ht into a long-standing problem of phase diagram and 



localization in f-electron systems 



nsignt 1 

HQ 



and has been used also to describe correlation effects 



2l\ . As an alternative to the above 



in half-metallic ferromagnetic materials like NiMnSb 
mentioned band structure methods, accurate self-consistent methods for solving the local 
Kohn-Sham equations based on LDA in terms of Green's functions have been developed 



y, y, Q- 



For that reason 



within the multiple scattering theory (KKR-method) |2z 
the KKR-method can be combined, as it will be shown below, in a natural way with the 
LDA+DMFT approach. A further appealing feature of this scheme is that the CPA alloy 
theory can also be incorporated very easily. 

In order to account within LDA-band structure calculations for correlations an improved 

fin 

hybrid Hamiltonian was proposed by Anisimov et al. 26J, 27]. In its most general form such 
a Hamiltonian is written as 



H = Elba + Hjj — H DC , (1) 

where Hlda stands for the ordinary LDA Hamiltonian, Hjj describes the effective 
electron-electron interaction and the one-particle Hamiltonian Hoc serves to eliminate dou- 
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ble counting of the interactions already accounted for by Hlda- 

Using second quantization a rather general expression for Hjj is given by: 



H u - 7, U ijkfiifilifinkC n l , (2) 



2 

n.ijkl 



where n runs over all the sites of the crystal R n and the creation (C^) and annihilation 
(C) operators are defined with respect to some subset of localized orbitals <pi{r — R n ). The 
constants U^ kl are matrix elements of the screened Coulomb interaction v(r — f '): 

U?W = j Rn)<f>}(r' - Rn)v(r- r')cj> k {r' - R n )Mr- R n )drdr' . (3) 



exactly, thus various meth- 

n 

^| . Among them one of the 



The resulting many-particle Hamiltonian can not be diagonalized_ 
ods were developed in the past to find an approximate solution 
most promising approaches is to solve Eq. (^J) within dynamical mean field theory, a method 
developed originally to deal with the Hubbard model. 

The main idea of DMFT is to map a periodic many-body problem onto an effective single- 
impurity problem that has to be solved self-consistently. For this purpose one describes the 
electronic properties of the system in terms of the one particle Green's function G(E), being 
the solution of the equation: 

(E - H - t(E))G = 1 , (4) 

where E is the complex energy and the effective self-energy operator S is assumed to be a 
single-site quantity for site n: 

±(E) = J2\<l>m)^(E)M. (5) 

ij 

Within DMFT, the self-energy matrix is a solution of the many-body problem of 

an impurity placed in an effective medium. This medium is described by the so called bath 
Green's function matrix Q defined as: 

gr j \E) = Gr j 1 (E) + ^(E), (6) 

where Gij(E) is calculated as a projection of G(E) onto the impurity site: 

G l3 {E) = (^GiE)^) . (7) 



As the self-energy depends on the bath Green's function Qij(E) the DMFT equa- 

tions have to be solved self-consistently. Accordingly, from a technical point of view the 
problem can be split into two parts. One is dealing with the solution of Eq. (0J) and the 
second one is the effective many-body problem to find the self-energy H^E). Within the 
present work, the first task is solved by the KKR band structure method, as described below 
in Sec. Ill Al The details of solving the many-body effective impurity problem based on the 
fluctuation exchange (FLEX) approximation will be presented in Sec. Ill Bl 

A. The KKR+DMFT formalism 

In this section we present an extension of the well known KKR equations in order to 
include the local, multi-orbital and energy dependent self-energy produced by the many- 
body solver (see section 111 Bj) . In the framework of the multiple scattering formalism the 
solution of Eq. (J3| is constructed in two steps. For the first step one has to solve the so 
called single-site scattering problem, to obtain the regular (Z) and irregular solution (J) of 
the corresponding Schrodinger (or in our case Lippmann-Schwinger) equations as well as a 
scattering amplitude expressed in terms of the single-site t-matrix. 

1. Solution of the single-site problem 

The solution of the single-site problem can be worked out easily in the same way as in 
the full-potential description jscj. This way one finally gets the single-site t-matrix for the 
LDA+DMFT case. In terms of the wave functions the single-site quasiparticle equation to 
be solved for each spin channel a reads 



In the following we omit the spin index a for the moment keeping in mind that for a spin- 
polarized system described in a non-relativistic way one has to solve Eq. (jHJ) for each spin 
channel independently. For the solution ty u (r) one can start from the ansatz: 






(9) 



L 
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where the partial waves ^lv{t) are chosen to have the same form as the linearly independent 
solutions for the spherically symmetric potential: 



*r„(r) = *Lv(r)Y L (r) , (10) 

with L = (I, mi) standing for the angular momentum and magnetic quantum numbers and 
Yi(f) are spherical harmonics. Inserting the ansatz (JSJ) into the single-site equation (jHJ) 
and integrating over angle variables leads to the following set of the coupled radial integro- 
differential equations: 

* Lv {r,E) = I r' 2 dr'Y, LL „{E) (f>i(r)<i>i»(/)*Lu(r',E) , 

L" J 

(11) 

For a general non-diagonal self-energy a similar radial equation (|11|) shall be written for 
the left-hand side equation. If one makes a rather natural choice of the localized subset of 
functions being just (phif) = (t>i{r)Yi,{f) ( see below). In principle these equations can be 
solved by summing a corresponding Born series. In this work, however, we simplified the 
equations taking advantage of the following special representation for the self-energy: 



/ d 3 r'i:(f,r\E)Mr',E) = J2 [ dV£^(£)4,(r^(f ')M? '> E ) « X) ^M E )Mr, E) 

J L L 

(12) 

This way the Eq. (fTTj) becomes a pure differential equation: 



d 2 1(1 + 1) 

dr 2 r 2 



V(r)+E 



V Lu (r,E) = Zll>(E) ^v v {r,E) . 



(13) 



After having solved the set of coupled equations for the wave functions one gets the 

n 

corresponding single-site t- matrix by introducing the auxiliary matrices a and b |3J: 



a Lu (E) = -ipr 2 [h L (pr),^ u L (r)} r 
b Lu (E) = -ipr 2 [hi(pr),^ L (r)] r . 



(14) 



Here p = \J~E~ is the momentum, h^(pr) are Hankel functions of the first and second kind 
and [. . ] r denotes the Wronskian. Evaluating the Wronskians at Wigner-Seitz radii rws one 



finally has |3j,|22|: 



t{E) = -{a{E)-b{E))b~\E). 



(15) 



The regular wave functions Z used to set up the electronic Green's function within the 
KKR-formalism j^J are obtained by a superposition of the wave functions ^ u according to 
the boundary conditions at r = r ws : 

Z L (r, E) = J2 Cl*„(r) ^ J>'(r, E)t(E)^ L , - ipht(f, E) , (16) 

v L> 

The irregular solutions Jl needed in addition are fixed by the boundary condition 

J L (f,E) r ^ s j L (f,E) (17) 

and are obtained just by inward integration with the functions being the spherical Bessel 
functions. 



2. The multiple scattering Green's function 

Having constructed a set of regular (Z) and irregular (J) solutions of the single-site 
problem together with the t-matrix the corresponding expression for the Green's function 
reads Q: 



G(r n + R n ,C + Rm,E) = J2 Z dr n ,E)T2 m L ,(E)Z*(C,E) 



L,L' 



- 6 nrn J2{M?n,E)JZ(CE)Q(K 

L 

+ J L (f n ,E)Zl(CE)Q(r n ~r;)}. 



(18) 



Here the superscript x is used to distinguish between the left and right hand solutions to 
Eq. (JBJ); i.e. for example \Z > and < Z x \ are solutions to the adjoint equations jsi^ : 



(H + t-E)\Z) = 
(Z x \(H + ± -E) = . 



(19) 
(20) 



The central quantity in Eq. (JTSj) is the scattering path operator r which for the case of a 
periodic crystal can be obtained from the Brillouin zone (BZ) integration: 



t2v(E) 



V, 



d 3 k 



BZ Jv BZ 



r 1 (E)-G(k,E) 



L'L 



(21) 
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where Vbz is the volume of the first Brillouin-zone and R nm = R n — R m with R n ( m ) denoting 
the lattice vector specifying the position of the unit cell n(m) and the matrix t~ 1 (E)—G(k, E) 
occurring in the integral is known as the KKR matrix. The matrix G(k, E) is the Fourier 
transform of the real space KKR structure constant matrix that depends only on the relative 
positions of scatterers. 

Given the local nature of the many-body solver used within the DMFT approach, the 
KKR Green's function (J18|) has to be projected accordingly to the matrix (see Eq. flj|)). 
The projection is performed through the following integration: 



Gl m v {E) = E (/ d^Hv^Z^E)^ t2?,(E) (I dV 2 Z* (r a ,E)0 L ,(r a )) 

SnmJ2( j d ^ (£ 1 dh^l^ZL^EYj Jl{T 2 ,E)4> L ,{Y 2 ) 

Li 

+ J d 3 r 2 QT" d'n^^J^E)^ Z^{v 2l E)<p LI {v 2 )) . (22) 

The impurity Green's function G^,(E) (actually G^jfiE) for both spin channels) repre- 
sents the input into the solution of the effective impurity problem presented below. As the 
DMFT-approach (see next section) concentrates on the correlation among electrons of the 
same angular momentum I only the I — /-subblock of this matrix will be used in the follow- 
ing. For the transition metal systems dealt here this implies that only the d — <i-subblock is 
considered with </>z,(r) being appropriate reference wave functions with / = 2. 



B. Solution of the effective impurity problem 



Our approach to achieve a solution of the many-body effective impurity problem is based 
on the fluctuation exchange (FLEX) approximation 28] but with a different treatment of 
particle-hole and particle-particle channels. The particle-particle channel is described by a 



T-matrix approach 



32 



33j giving a renormalization of the effective interaction, the latter 



one being used explicitly in the particle-hole channel 



16 



34] 



The symmetrization of the bare U matrix is done over particle-hole and particle-particle 
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channels: 



j id — UP — TP 

jjm = 

U s = -(W +U i ) 

U l = -(IP -U 1 ) . 

As indicated above, here and in the following only matrix elements with respect to the rf-like 
reference wave functions 4>l have to be considered. The above expressions are the matrix 
elements of bare interaction which can be obtained with the help of the pairwise operators 
corresponding to different channels: 



particle-hole density 



• particle-hole magnetic 



particle-particle singlet 



particle-particle triplet 



l 12 



-^=(4^ + 0+^21) 



(23) 



m+ = c+ T c 2i (24) 
m i2 = c u c 2T 



S12 = -^=(C U C 2T -C 1T C 2i ) 

S12 = —^(c^c^-c^c^) (25) 



*12 = ^|( C H C 2T + C 1T C 2|) 

t 12 = — 7=(cjpcjj_ + c^c^) 

= C 1T,I C 2|,T 

*i2 = c ir,4 c 2i,T ■ (^6) 



10 



These operators describe the correlated movement of the electrons and holes below and 
above the Fermi level and play an important role in defining the spin-dependent effective 
potentials W™ m2msm4 . The one-electron Green's function matrix containing the many-body 
interaction, described by the self-energy 'E mm ' a (iu n ) is given by the Dyson equation 

Gmln 'ai^n) = (^n + f)^mm ' ~ h mm > a — S mm > a (iuJ n ) (27) 

where /i is the chemical potential, u n = (2n + 1)tt/(3 are Matsubara frequencies and /3 = 1/T 
is the inverse temperature. The GW type of diagrams are summed up self-consistently to 
produce the self-energy. For getting the self-energy we use a two-step FLEX approximation. 
This means that first of all the bare matrix vertex is replaced by the T-matrix approach 



32, 



which will be used in the calculation of the particle-hole channel. In the Kanamori 
T-matrix approach the sum over the ladder graphs may be carried out with the aid of the 
so called T-matrix which obeys the Dyson-like integral equation: 

< 13|T ,J<T '(^)|24 > = 

<13M24> - ^^2^2 < 13\v\57 > g^{iu)g^ {iQ - iu) <68\T aa '{iQ)\2A> . 

P u 5678 

The Hartree and Fock contribution are obtained replacing the bare interaction by a T-matrix: 

= ^^<13|T^'(^)|24>^(^-H (28) 

" n 34ct' 

= -7? EE < 141^(^)132 > %(<0-iw) . (29) 

" Q 34 

In the low-density limit the self-energy should be the summation over diagrams for repulsion 
of two holes below Ep (ladder approximation). Going beyond the low density limit means 
the inclusion of excitations of electrons from states below the Fermi level into the unoccupied 
part of the d band. This process renormalizes the hole states below Ep and put new poles 
for the Green's function. 

Combining the density and the magnetic parts of the particle-hole channel we can write 
the expression for the interaction part of the Hamiltonian [nJI3 : 

H v = ^Tr{D + * V 11 * D + m + * * m~ + rrT * * m + ) , (30) 

where D is a row matrix with elements (d,m°), and D + is a column matrix with elements 
(d + rriQ ). We denote by * matrix multiplication with respect to the pairs of orbital indices. 

11 



The expression for the effective potential is: 

1 / T/dd T/dm \ 

V(i») = \ / (3D 
(V^)i 2 34 = <13|T^|42> . (32) 
The matrix elements of the effective interaction for z or longitudinal spin-fluctuations are: 
V dd = < 13 I T<J<T 'I 42 > - < 13|T a ' CT '|42 >) 

2 a a' 

V dm = V md = l^2^(< 13|T CT<T |42 > - < f3|T aCT |24 > + < I3|T CT ' a |42 >) 

aa ' 

V mm = ^(^^'<i3|T™'|42>-<I3|T ctV '|42>). 

a a 1 

For finite temperature the definition for the spin dependent Green's function is: 

QUr) = - <T T c la (r)ct(0) > 

Jo 

The corresponding expressions for the generalized longitudinal x" an d transversal x 1 " sus- 
ceptibilities are: 

X \iu) = [1 + V^T n (iuj)Y l * rK(iu) (33) 
= [1 + V 11 * XoM]" 1 * xlM , (34) 

where T(iu) represent the Fourier transform of the empty loop: 

^mira2m3m4 f''") Qmimz (^~)^m4mi ( (35) 

and the matrix of the bare longitudinal susceptibility is: 

i / r TT + r u r TT - r u \ 

x » M = 2^ r „_ rU m+ruj ' (36) 

The four matrix elements of the bare longitudinal susceptibility represent the density- 
density (dd), density-magnetic (dm ), magnetic-density (m°d) and magnetic-magnetic chan- 
nels (m°m°). The matrix elements couple longitudinal magnetic fluctuations with density 
magnetic fluctuations. In this case the particle hole contribution to the self-energy is: 

SSW^^W^'W, (37) 



34a > 
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with the particle-hole fluctuation potential matrix 



and the spin-dependent effective potentials defined as: 



W aa '(iu) =1 TT 11 I , (38) 





2 


* (x 11 - 


Xo) * 


vw 




w u = 


2 


* (x 11 - 


xS) * 






w n = 


V 

2 m 


*(x + ~ 


-Xo 


-)* 




w u = 


2 m 


*(x" + 


- Xo 


+)* 





(39) 

The definitions for jfi and Xo differ from those of x" and xJl ) respectively, by the replacement 
r^T <-> in Eq. ( )3f)jl . The complete expression for the self-energy is finally given by: 

S = £< TH ) + £< TF ) + E« . (40) 

The attractive feature of the present approach is that it leads to an exact expression 
for the self-energy in the limit of a small number of holes in the d band. These conditions 
are satisfied with high accurac y in the case of Ni. Further details and justifications of this 
approach can be found in Ref. |3 



C. Treatment of disordered alloys 

In this section we review the KKR-CPA approach and present a simple and transpar- 
ent electronic theory that combines the treatment of disorder and correlation on the same 
footing. After several decades of intense research the problem of interacting electrons in dis- 
ordered alloys still induce numerous investigations both experimentally and theoretically. In 
the weakly disordered limit [35[ both disorder and interaction can be treated in a perturba- 
tive way; note that this perturbation theory is not trivial, in particular, a non-Fermi-liquid 
behaviour appears. For strong disorder Anderson localization effects eventually lead to the 
breakdown of the metallic phase and a metal-to-insulator transition takes place (for a re- 
. It was realized recently that the Hubbard model can be solved exactly 
in the limit of infinite space dimensionality d = 00 and in this case the Mott metal-insulator 
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transition can be described in the framework of dynamical mean-field theory (for a review, 
see Ref . U) . The presence of disorder in d = oo increases the complexity of the problem: the 
cavity field varies from site to site reflecting the random environments in which a given site 
is embedded JstJ. Fortunately, for d = oo the problem can be simplified due to a (infinitely) 
large number of neighbours, in this case the cavity fields become independent of disorder 
and only local disorder fluctuations survive. We will adopt this approach which is flexible 
enough to allow for study numerous interesting questions in connection with an interplay 
between correlations and local disorder Furthermore it is supported by the arguments 
given by Drchal et al. [5]. These authors pointed out that an averaged coherent poten- 
tial for disordered interacting systems can be constructed using the so-called terminal-point 
approximation. Using a local mean-field approximation to treat electron correlations, the 
corresponding self-energy gets diagonal in the site representation. This allows to use the 
coherent potential alloy theory (CPA) j3] for the configurational averaging in the usual 
way. 

Among the electronic structure theories, those based on the multiple scattering formalism 
are the most suitable to deal with disordered alloys within the coherent potential approxima- 
tion (CPA). CPA is considered to be the best theory among the so-called single-site (local) 
alloy theories that assume complete random disorder and ignore short-range order (^J • Com- 
bining the CPA with multiple scattering theory leads to the KKR-CPA scheme, which is 
applied nowadays extensively for quantitative investigations of the electronic structure and 
properties of disordered alloys j^lS]- Within the CPA the configurationally averaged prop- 
erties of a disordered alloy are represented by a hypothetical ordered CPA-medium, which 
in turn may be described by a corresponding site-diagonal (n = m) scattering path operator 
t cpa rj-^g corresponding single-site t-matrix t CPA and multiple scattering path operator 
r CPA are determined by the so called CPA-condition: 

x A r A + x B r B = r CPA . (41) 

Here a binary system A^B!^ composed of components A and B with relative concentrations 
xa and xb is considered. The above equation represents the requirement that embedding 
substitutional^ an atom (of type A or B) into the CPA medium should not cause additional 
scattering. The scattering properties of an A atom embedded in the CPA medium, are 
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represented by the site-diagonal (n = m) component-projected scattering path operator r A 

r A = r GPA [l+te 1 -tSA)r CPA ]- 1 , (42) 

where t& and ^cpa are the single-site matrices of the A component and of the CPA effective 
medium. A corresponding equation holds also for the B component in the CPA medium. 
The coupled sets of equations for r CPA and t CPA have to be solved iteratively within the 
CPA cycle. 

It is obvious that the above scheme can straightforwardly be extended to include the 
many-body correlation effects for disordered alloys. As was pointed out in Sec. Ill Al within 
the KKR+DMFT approach the local multi-orbital and energy dependent self-energy (Ea{E) 
and Tjb{E)) is directly included in the single-site matrices and tg, respectively. Having 
solved the CPA equations self-consistently, one has to project the CPA Green's function 
onto the components A and B by using Eqs. (|22j) and (J42j) . In Eq. (J22j) the multiple 
scattering path operator t^ l ,(E) has to be replaced by the component-projected scattering 
path operator r A £7 of an A-atom in a CPA medium. The components Green's functions 
Gi=A,B are used to construct the corresponding bath Green's functions for which the DMFT 
self-consistency condition is used according to Eq. (fU): 

G; = \ B (E) = G;} AB (E) + Z i=AB (E) . (43) 

The many-body solver presented in section 111 ±31 in turn is used to produce the component 
specific self-energies T li= A,B{.E): 

Vi=A, B {E) = ^i=a,b(E) [Qi=A,B{E)\ . (44) 
D. The self-consistency cycle 

Finally a description of the flow diagram of the self-consistent LDA+DMFT approach 
is presented in Fig. ((TJ). The radial equation Eq. [TT] provides the set of regular (Z) and 
(J) irregular solutions of the single-site problem. Together with the t matrix, the scattering 
path operator r Eq. (J21|) and the KKR Green's function is constructed Eq. To solve the 
many-body problem the a projected impurity Green's function is constructed according to 
Eq. [221 The LDA Green's function G^,(E) is calculated on the complex contour which 
encloses the valence band one-electron energy poles. The Pade analytical continuation is 
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used to map the complex local Green's function G2l' (-^0 on the Matsubara axis which is used 
when dealing with the many-body problem. In the current implementation the perturbative 
SPTF (spin-polarized T-matrix + FLEX) solver of the DMFT problem described above 
is used. In fact any DMFT solver could be included which supplies the self-energy £(u;) 
as a solution of the many-body problem. The Pade analytical continuation is used once 
more to map back the self-energy from the Matsubara axis to the complex plane, where 
the new local Green's function is calculated. As was described in the previous sections, 
the key role is played by the scattering path operator r£ which allows us to calculate 

the charge at each SCF iteration and the new potentials that are used to generate the new 
LDA Green's function. In practice it turns out that the self-energy converges faster than 
the charge density. Of course double counting corrections have to be considered explicitly 
when calculating the total energy (not done here). Concerning the self-energy used here the 
double counting corrections are included when solving the many-body problem (see Ref.Q. 



III. RESULTS AND DISCUSSION 



To demonstrate the capability of our approach we first applied it to the 3d metals Ni 
and Fe. Although these metals are more or less adequately described in the framework of 
standard LDA, nevertheless, there are some features in the experimental properties which 
are due to correlation effects that are not adequately described on this basis. In addition, 
there are numerous investigations in the literature that seek for an improved description of 
correlation effects in these systems and that can be compared with. 



A. Numerical details 



The self-consistent LDA+DMFT calculations were carried out for the experimental 
ground state crystal structures, i.e. fee for Ni, bec for Fe and fee for Ni rich Fe^Mi^ 
alloys. The lattice parameters were fixed at the experimental values (Fe: 5.406 a.u., Ni: 
6.658 a.u., for Fea;Nii_ x .: see RefQ). The Green's function was calculated for 32 energy 
points distributed over semicircular contour. The Brillouin zone integration has been per- 
formed on a uniform grid, taking into account the symmetry of the system. As a suitable 
reference wave functions 0l(^* — R) we have choosen a radial solution of the Schrodinger 
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equation for the spherically symmetric LDA non-magnetic potential, that is determined for 
an appropriate energy (E=0.7 Ry). The DMFT parameters, average Coulomb interaction 
U, exchange energy J, and temperature T used in the calculations are listed in Table 



B. Results for bcc-Fe and fcc-Ni 



To demonstrate the applicability of the scheme presented above band-structure calcula- 
tions for bcc-Fe and fcc-Ni have been performed. The results of LDA+DMFT calculation for 
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both systems have been already several times discussed in detail in the literature 

The density of states curves resulting from a plain LDA and a LDA+DMFT calculations 
are shown in Figs. El and H3 for Fe and Ni, respectively. For the LDA+DMFT calculations 
we used the DMFT parameters as given in Table [IJ The density of states curves for Fe 

n 

and Ni are in reasonable agreement with corresponding previous LMTO+DMFT [5J, as well 
as EMTO+DMFT |l2j calculations. The same is true also for the spin magnetic moments 
(see Table H}. The spin magnetic moments are some what higher in comparison with the 
EMTO+DMFT results, jl^] From Figs 121 and El one can see that in bcc-Fe the correlation 
effects are much less pronounced than in fcc-Ni. This is due to the large exchange splitting 
for Fe and the bcc-structure dip in the minority density of states [4lJ. In the case of Ni 
the LDA+DMFT calculations account for all expected influences of the density of states in 
satisfying way. As can be seen from Fig. El the density of states reflects all three main corre- 
lation effects: the 30% narrowing of the occupied part of the d-band, about 40% decrease of 
exchange splitting and the presence of the famous 6eV satellite compared to the LDA DOS. 
However, the position of the 6eV sattelite is shifted somewhat to lower binding energies. 
This shift and the large broadening of the resonance is due to the perturbation approach of 
the DMFT solver of the effective impurity problem used here, 



C. Results for fcc-Fe x Nii_ x . disordered alloy 

As mentioned above, the scheme presented here allows in a straight forward way to deal 
with disordered alloys. To demonstrate how this works we carried out a set of LDA as well 
as LDA+DMFT calculations for fcc-Fe^Nii.^ disordered alloy for various concentrations. 
For the LDA+DMFT we used the same DMFT parameters (U, J and T) as in the case of 
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pure bcc-Fe and fcc-Ni (see Table UJ). 111 Fig. H] the element resolved as well as the total spin 
magnetic moments are shown. Although the difference between LDA and LDA+DMFT 
moments for Fe is rather small one can see an interesting trend. In contrast to the pure 
bcc-Fe case the LDA+DMFT moments for Fe in fcc-Fe^Nix^ alloy are slightly larger than 
corresponding LDA ones. In the case of Ni, on the other hand, a decrease of the magnetic 
moment was obtained as in the case of pure fcc-Ni (see Table UJ). Comparing the average 
moments in Fig. 0] with experiment one finds rather good agreement already for LDA- 
based spin moment. In spite of its present limitations, the LDA+DMFT scheme does not 
spoil the overall behaviour for the concentration dependence of magnetic moments. 

Finally, in Fig. the resulting concentration dependence of the self-energy is shown. 
We present the results for the real part of the self-energies for Fe and Ni atoms for the 
t2 9 symmetry (the results for the e g symmetry are similar and hence not plotted). It is 
interesting to note that the slope of the self-energy near Fermi level Z = j^\ e=Ef which 
defines the mass renormalisation and leads to the narrowing of the band practically does not 
depend on the concentration. On the other hand, for the high energy part of the self-energy 
one sees rather noticeable differences giving raise to the changes in the sattelite structure. 

IV. SUMMARY 

A scheme has been presented that allows to combine the KKR band structure method 
and the LDA+DMFT approach to deal with correlated systems. Its applicability has been 
demonstrated by results for ferromagnetic bcc-Fe and fcc-Ni. For this systems a good 
agreement with previous LDA+DMFT methods has been found. In addition we combined 
LDA+DMFT scheme with the CPA to deal with disordered alloy. As an example we pre- 
sented results for fcc-Fe^Nix^ disordered alloy. 
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EF 



FIG. 1: The complex energy contours used within the selfconsistent LDA+DMFT approach, com- 
bined with the KKR formalism. The solution of the radial equation allows the evalaution of the 
single site scattering matrix t(E) and the scattering path operator t^jj{E) from which the KKR 
Green function is constructed. The projection of the Green function is perfomed according to 
Eq. (|'22[). The impurity Green function is then used to solve the many-body problem within the 
spin polarized T-matrix FLEX solver of the DMFT approach. 
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FIG. 2: Spin resolved density of states of bcc-Fe as calculated within LDA(dashed line) and 
LDA+DMFT (full line) using the KKR-method. (DMFT parameters: U=2eV, J=0.9eV, T=400K) 



TABLE I: The DMFT parameters average Coulomb interaction U, exchange energy J and tem- 
perature T used in the calculations for bcc-Fe, fcc-Ni and fcc-FesoNiso- In addition the theoretical 
spin magnetic moments as calculated by the LDA and the LDA+DMFT methods are shown for 
bcc-Fe and fcc-Ni. Magnetic moments for fcc-Fe x Nii_ x alloy are presented in Fig. 0] 
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FIG. 3: Spin resolved density of states of fcc-Ni as calculated in the LDA (dashed line) and 
LDA+DMFT (full line) using the KKR-method. (DMFT parameters: U=3eV, J=0.9eV T=400K) 
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Fe content (at. %) 



FIG. 4: Spin magnetic moments of Fe (circles) and Ni (squares) in Fe x Nii_ x alloy calculated using 
plain LDA (open symbols) and the LDA+DMFT method (full symbols). 
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FIG. 5: Left: Concentration dependence of the real part of the spin resolved self-energy for Fe. 
Only results for t2 g d-orbitals are shown. Right: Same as in the left panel but for Ni 
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